TECHNISCHE UNIVERSITAT MUNCHEN 



Dynamic scattering from solutions of semifiexible polymers 

Klaus Kroy and Erwin Frey 
Institut fiir Theoretische Physik, Physik-Department der Technischen Universitdt Miinchen, 
James-Franck-Strafie, D-85747 Garching, Germany 
(submitted to PRE) 

The dynamic structure factor of semifiexible polymers in solution is derived from the wormlike 
chain model. Special attention is paid to the rigid constraint of an inextensible contour and to 
the hydrodynamic interactions. For the cases of dilute and semidilute solutions exact expressions 
for the initial slope are obtained. When the hydrodynamic interaction is treated on the level of a 
renormalized friction coefficient, the decay of the structure factor due to the structural relaxation 
obeys a stretched exponential law in agreement with experiments on actin. We show how the 
characteristic parameters of the system (the persistence length Ip, the lateral diameter a of the 
molecules, and the mesh size of the network) are readily determined by a single scattering 
experiment with scattering wavelength A obeying a <^ \ <^ (.p and \ < £,m- We also find an 
exact explicit expression for the effective (wave- vector-dependent) dynamic exponent z{k) < 3 for 
semifiexible polymers and thus an enlightening explanation for a longstanding puzzle in polymer 
physics. 
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I. INTRODUCTION 

Recently, there has been increasing interest in biologi- 
cal materials research [1] . The physical properties of col- 
loids, liquid crystals, and macromolecular networks are 
of prime importance for the structure and function of bi- 
ological entities such as cells and muscles. On the other 
hand, biology provides physicists with some of the most 
pertinent model systems to test their theories of soft mat- 
ter. Among these systems we will concentrate on solu- 
tions of semifiexible macromolecules here. The challeng- 
ing problems associated with semiflexibility, which have 
been attacked by numerous groups over many years, have 
received recent attention. Semiflexibility has been recog- 
nized as a crucial property for the understanding of many 
peculiar features of DNA [2,3] and actin [4-8]. This may 
not be too surprising in the case of actin, which poly- 
merizes into filaments that are rarely much longer than 
their persistence length. It is perhaps less obvious for 
DNA molecules, which are typically orders of magnitude 
longer than their persistence length and thus have an 
overall flexible appearance. 

In this paper we will concentrate on the dynamic as- 
pects of semiflexibility. In contrast to the dynamics of 
flexible polymers [9] the dynamic properties of semifiexi- 
ble polymers are still not very well understood: Whereas 
a simple, analytically tractable basic model for flexible 
polymers has been known for a long time, we lack such a 
generally accepted simple model in the case of semifiex- 
ible polymers. It seems worth mentioning the reason for 
simplicity and complexity in both cases. The standard 
model for the statics and dynamics of a flexible polymer 
is the Gaussian chain [9], which represents the connec- 
tion of the monomers by an isotropic harmonic poten- 
tial characterized by a single parameter, the mean square 
end-to-end distance of the polymer. This description re- 



produces many of the universal large scale properties of 
flexible polymers, which are purely entropic in origin. 
The universality may be understood as a consequence of 
the central limit theorem. The harmonic theory allows 
extensions to nonideal problems by perturbation theory. 
The simplicity of the model is due to the fractal struc- 
ture of the Gaussian chain: it looks exactly the same on 
all scales. This property is clearly not shared by real 
polymers. They are coiled on large scales, but they are 
rodlike on length scales below their persistence length ip. 
A more realistic polymer model than the Gaussian chain 
is the so-called "Kratky-Porod" model [10]. In this model 
the conformation is derived from an effective free energy 

111], 

which takes into account the energy cost of bending the 
contour. This is the contour integral over the square of 
the local curvature multiplied by the bending modulus 
K. The conformation resulting from the above free en- 
ergy together with the rigid constraint \drs/ds\ = 1 of 
an inextensible contour is known as the "wormlike chain". 
It is not a fractal as for the Gaussian chain model but 
a differentiable curve, which is indeed rodlike on short 
distances and coiled on large scales. Short and long dis- 
tances are measured with respect to the decay length of 
the tangent-tangent correlations, the persistence length 
£p = t^/keT [12]. 

If one does not look for a description of the specific mi- 
croscopic details of a macromolecule but wants to under- 
stand global universal properties, shared by large classes 
of molecules, then one need not worry about the micro- 
scopic discrepancy between a real polymer and the Gaus- 
sian chain model in the case of proper flexible polymers 
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with ^p of the order of the lateral diameter a. In this ease 
the local rodlike structure is part of the nonuniversal mi- 
croscopic details. However, £p/a is somewhat larger than 
1 in several important cases, e.g., for polystyrene, DNA, 
actin etc. (for actin tp/a ~ 10^). One may still hope 
that the intrinsic stiffness is negligible, when the total 
contour length L of the molecules is much larger than (.p, 
because the large scale properties are then dominated by 
the coil structure. 

As far as the dynamics is concerned, this heuristic ar- 
gument has one flaw with possibly serious consequences: 
hydrodynamic interactions play a crucial role for the dy- 
namic properties of polymers in solution, as they do for 
any system of Brownian particles in hydrodynamic sol- 
vents. Because the dynamics of the solvent, which me- 
diates these interactions, is usually much faster than the 
Brownian dynamics, the hydrodynamic interactions in 
dilute and semidilute solutions can be subsumed into an 
(instantaneous) mobility matrix, 

called the Oseen tensor, which is calculated from the 
Navier-Stokes equation for the pure solvent [9]. Multi- 
plied by a force acting at the origin it gives a contribution 
to the velocity v of the solvent at any point in space r. 
This velocity field has the characteristic feature to decay 
like 1/r in real space, where r is the length of the vector 
r; i.e., it is both long ranged and singular at the origin. 
Physically speaking, the singularity is a consequence of 
idealizing a physical particle as a mass point and has to 
be cut off at about the diameter a of the polymer. For 
a real polymer, this short distance divergence has the 
consequence that its local semiflexible structure, which 
one would possibly prefer to neglect, may markedly pro- 
nounce itself in the dynamic properties, even if the poly- 
mer as a whole looks rather flexible. If one decides to 
take care of this effect, one immediately runs into a prob- 
lem: The longitudinal degree of freedom of any contour 
element of the polymer is suppressed by the presence of 
its neighboring contour elements. This rigid constraint 
of a locally inextensible contour is the source of the dif- 
ficulty in modeling the dynamics of (intrinsically) semi- 
fiexible polymers. It renders awkward any general theory 
[11,13,14] which tries to represent this property faithfully. 
On the other hand, models that relax the constraint too 
much — as e.g., the so-called Harris-Hearst-Beals model 
[15] and its latest descendants [16 19] - include artificial 
stretching modes and find a Gaussian distribution for all 
spatial distances along the contour; i.e., the essence of 
semiflexibility has obviously been lost. The correct ra- 
dial distribution function of a semiflexible polymer with 
L K £p is actually very different from a Gaussian dis- 
tribution [20]. It is not peaked at the origin and hence 
cannot be approximated by a Gaussian form. Moreover, 
these "Gaussian" models treat the thermal fluctuations 
to be isotropic. However, as was pointed out by many 



workers in the field (see, e.g., Refs. [14,21]), it is concep- 
tually important to be aware of the local anisotropy of 
the bending undulations caused by the rigid constraint 
of constant contour length. 

There is a special ease, where semiflexible dynamics 
can be treated analytically with moderate expense. This 
is the limit of a weakly bending rod, which also has been 
addressed previously [21-23]. However, these workers 
did not arrive at the results given below. Especially the 
crucial point of the local anisotropy of the undulations, 
which was recognized in Ref. [21] but not discussed in 
Ref. [23], will be analyzed more closely in this contri- 
bution. We will show that this complication does not 
necessarily prevent an analytical approach. In a forth- 
coming publication [24] we will also give a quantitative 
comparison of our results for the decay of the structure 
factor with new experimental data. 

Below, we will demonstrate how one can incorporate 
the rather complex hydrodynamic interactions and ef- 
fects from chemical cross-linking in a semidilute solution 
into the description on a reasonable level of accuracy. To 
be specific, for the computation of the dynamic structure 
factor we require the scale separation 

a <^ X<^ip, L and X<^rn (1-3) 

to be realized. Here we have introduced the symbol 
for the mesh size in a semidilute solution and the scatter- 
ing wavelength A = 27r/fc. The latter should be shorter 
than the mesh size to probe the dynamics of single poly- 
mers. The conditions in Eq. (1.3) are generic for such 
important cases as neutron scattering or light scattering 
experiments from semidilute solutions of DNA or actin, 
respectively. 

The outline of the paper is as follows. In Sec. II 
we derive the functional form of the decay of the dy- 
namic structure factor for semiflexible polymers in so- 
lution. Special attention has to be paid to the rigid 
constraint of constant contour length and to the hydro- 
dynamic interaction. We demonstrate how this can be 
achieved in the case of interest. We give the exact an- 
alytical result for the time decay of the structure factor 
and show that it can be reduced to a simple stretched 
exponential law, which is readily applied to determine 
various parameters of interest from experimental data. 
In Sec. Ill we calculate the initial slope of the structure 
factor. This can be used to measure the microscopic lat- 
eral diameter a of the polymer by a dynamic scattering 
experiment with A o. We show that the short time 
dynamics of semiflexible polymers is characterized by a 
wavelength-dependent dynamic exponent z(k), which is 
computed explicitly. What is remarkable about our re- 
sults compared to other predictions available in the litera- 
ture [25,21] is that they are simple analytical expressions, 
which at the same time excellently fit experimental data 
[6]. Finally, in Sec. IV we summarize our findings and 
outline some possible experimental applications. 
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II. STRUCTURAL RELAXATION 

For the following wc suppose Eq. (1.3) holds. We con- 
centrate on the common case of a semidilute solution, the 
dilute limit being included as a special case. With the 
conditions in Eq. (1.3) the decay of the structure factor is 
mainly caused by the local conformational fluctuations, 
whereas the global structure of the polymer network stays 
virtually fixed during the characteristic decay time. This 
allows for a description of the dynamics purely in terms 
of the transverse undulations of a weakly bending rod, 
i.e. in terms of the Langevin equation 
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for the polymer contour r^. The random force (per 
length) is assumed to represent Gaussian white noise. 
The conformational energy is given by Eq. (1.1). The 
mobility matrix 
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serves to mediate the hydrodynamic interaction and to 
project out the forbidden longitudinal motion of the seg- 
ments. This is a convenient method to enforce the rigid 
constraint of fixed contour length mentioned above on 
the local dynamics. Please note that Eq. (2.2) there- 
fore differs from the Oseen tensor given in Eq. (1.2) by 
the opposite sign of the projector |'r)(r|. Whereas (as 
a consequence of the incompressibility of the solvent) 
the longitudinal direction is weighted twice in the usual 
Oseen tensor, it is now completely suppressed. (Re- 
call that we can neglect the center of mass and rota- 
tional modes as a consequence of the time scale separa- 
tion induced by Eq. (1.3).) The exponential prefactor 
in Eq. (2.2) accounts for the screening of the hydrody- 
namic self-interaction of a single molecule by the sur- 
rounding network: hydrodynamic interactions over dis- 
tances longer than the mesh size are largely suppressed 
because of the larger effective viscosity of the solution 
as compared to the pure solvent. A formal derivation of 
the exponential screening term has been achieved before 
within a self-consistent approach known as the "effective 
medium theory" [9] . One expects the screening length to 
be roughly equal to the mesh size. It was indeed shown 
to obey basically the same scaling law as a function of 
the concentration [26] . Simple geometrical considerations 
[27] suggest that the mesh size of a semidilute solution of 
semiflexible polymers with a large ratio tp/a obeys to a 
very good approximation the scaling law 



Cm C ^ . 



(2.3) 



until the mesh size becomes considerably larger than the 
persistence length. 



For the following wc will neglect the nonlocal nature 
of the mobility matrix Eq. (2.2) and replace it by the 
inverse of the effective transverse friction coefficient (per 
length) 



(2.4) 



for the transverse undulations. This coefficient may be 
obtained from Eq. (2.2) by taking the terms in parenthe- 
ses in Eq. (2.1) out of the integral and averaging over all 
segments s for a virtually straight contour. However, we 
want to point out briefly how one can proceed more rigor- 
ously. Working with the full mobility matrix Eq. (2.2) in 
the equation of motion, Eq. (2.1), one needs its Fourier 
transform. For the transformation to make sense, the 
short distance singularity mentioned above has to be cut 
off at the diameter of the molecule a. For the transverse 
undulations we thus set the transverse mobility equal to 
the prefactor in Eq. (2.2) with a short distance cutoff 
represented by another exponential screening term: 



fl±(r) = 



(2.5) 



This leads to a convenient regularization for the familiar 
[28,23,25] logarithmic mode number dependence of the 

hydrodynamic interaction: The Fourier transform of H± 
is diagonal with the diagonal elements H±{q) given by 



87r^fil(,) = -ln(-^ ^ -In 



For modes of wavelength longer than the screening length 
{Q^h ^ 1) the last formula reduces to the mode mobil- 
ity (2^x)~^ derived from Eq. (2.6) by the normal mode 
analysis (see the Appendix), but it predicts a somewhat 
smaller mobility H±{q) — Inqa/Snrj + 0{{q^h)~^) for 
the short wavelength modes with q^h ^1- It is not a 
critical approximation to neglect the weak wavelength de- 
pendence of the mobility in the following considerations. 
This will allow for a simple expression for the decay of 
the structure factor, and at the same time it captures 
the main effect of the hydrodynamic self-interaction of 
the single polymers as well as their mutual interaction. 
The error made in this approximation can in part be com- 
pensated by renormalizing to an "effective hydrody- 
namic screening length" . Since the original was a phe- 
nomenological parameter, which had to be determined 
from experiment anyway, this is not a serious shortcom- 
ing from a practical point of view. For the rather dilute 
case {k^m ^ 1) a theoretical estimate for the param- 
eter in Eq. (2.4) — motivated by the exact result, 
Eq. (3.4), obtained in Sec. Ill for the initial decay rate 
— is k^h — e^^® (see the Appendix). A further refinement 
would be tedious and its effects presumably undetectable 
in typical scattering experiments. 
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With the above considerations in mind we may use a 
"Rouse-like" linear equation with a renormalized friction 
coefRcient (± for the transverse local undulations rj-{t): 



(2.6) 



The transverse random force (per length) represents 
Gaussian white noise. Eq. (2.6) can be solved by a nor- 
mal mode analysis. Physically, from Eq. (1.3) one expects 
that the problem at hand should not depend sensitively 
on the special choice of the boundary conditions. This 
may be justified by a formal argument as well. The reader 
may generally (but not invariably) think of the normal 
modes as cosine functions (see the Appendix). Effects 
from chemical cross-linking of the polymers (clamped 
ends) will be incorporated later on. 

To calculate the decay of the dynamic structure factor, 
we start from the definition 

*) = ^ ^(exp[ife- (r„(t) - r„(0))]) (2.7) 

for the dynamic structure factor of a statistical system 
composed of N equal scattering centers (monomers). As 
we have explained above, our main interest is in situa- 
tions where the decay of the structure factor is mainly 
caused by the fast local bending undulations while the 
global structure of the network stays virtually fixed. 
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Therefore, we perform the conformational average in two 
steps: First we take the thermal average (. . .)t over the 
transverse undulations of a weakly bending rod keeping 
the mean orientation of the rod fixed in space. Then 
we average over an isotropic ensemble of rod orienta- 
tions (. . .)o- As the transverse undulations obey a linear 
Langevin equation, they have Gaussian correlations as 
the random forces. Using the abbreviations 

:= a{n - m) , <^(i) := {{r^{t) - ri{0)f)^ 

we can rewrite g{k, t) as 

'ik'--{r^{t)-ri{0)) 



-E 



exp 



exp [I 



and, after performing the thermal average, as 
1 ^ (exp {-k^'^niKt)!^] exp[iA;ll7eL] 



(2.8) 



Performing now also the oricntational average, we obtain 
the explicit general expression 



(2.9) 



with 



(2.10) 



To proceed further we introduce the Rouse-like decay 
time of the mode of wavelength 2L 

.. = ^(^)'. (2.11) 

which is immediately read off by dimensional analysis 
from Eq. (2.6) as the characteristic time scale. As an- 
other natural abbreviation we introduce the decay rate 



(2.12) 



which will be shown to govern the time decay of the struc- 
ture factor. The decay rate niay also be predicted by a 
simple scaling argument [28,23]: For the wormlike chain 
the amplitudes of small transverse undulations scale as 



„-L2 



oc r" /£p with the wavelength r''. Substituting the 
scattering wavelength A for the amplitude r-*- of the un- 
dulations and r'l for the contour length parameter L in 
Eq. (2.11) one obtains Eq. (2.12). This heuristic argu- 
ment is a special case of a more general theorem derived 
in Ref. [28] for the dispersion relation of the decay rate: 



7fc 



A:2/° 



with 



2C 



D + 2( + a 



the stretching exponent. Here ( denotes the roughness 

exponent (C = 3/2 for the weakly bending rod), D is the 
dimension of the fluctuating manifold, and a character- 
izes the hydrodynamic interaction (a = in the present 
case). 

After solving Eq. (2.6) by a normal mode analysis and 
some technical manipulations (see the Appendix) we find 



k^K^it) = -{lkt)il{xm,y) + k^KliO), (2.13) 



where we have introduced for t > the dimcnsionless 
variables ■= (t/TL^Y^'^, y := kTllm/[{lkt)Hk£p)i] 
and [29] 

/(.„,,):= p. ^(l-e-^^). (2.14) 



The lower limit of the integral serves to exclude modes 
of wavelength longer than 2Lm- By Lm we denote the 
contour length between adjacent clamped points along 
the contour. This infrared cutoff is a heuristic way to 
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take into account the effects of steric constraints, so- 
called entanglements. Strictly speaking, these entangle- 
ments cannot fix the polymers in space efficiently at short 
times (7^^^(t) ^ ^^), i.e., until an average contour el- 
ement has moved a distance of the order of the mesh 
size ^m- Only for times longer than 1/7^-1 do the con- 
straints suppress modes of wavelength longer than the 
entanglement length L„i = (3/2^^ ^p)-'^/^ [30]. However, 
this subtle point can safely be neglected, because Xm «C 1 
for short times (7^-1 1 <^ 1) anyway. If kL^ is sufficiently 
large the modes wTiich contribute most to the decay are 
not substantially disturbed until the structure factor has 
completely decayed, i.e., ^uTl^ — (k^m)^^^ ^ 1 and we 
can take the limit Xm — > in Eq. (2.14). On the other 
hand, if the condition k^m ^ 1 is not fulfilled, the steric 
constraints become important and the structure factor 
decays in two steps: For short times (7^^^(t) ^^) the 
decay is mainly due to the fast single chain dynamics. 
As a consequence of the infrared cutoff the struc- 
ture factor does not decay to zero but only to a value 
g{k,t > TLj/g{k,0) « exp[-(7feTi„)3/4/37r]. The fur- 
ther decay is due to the slower collective modes of the net- 
work, which are not considered here. An obvious question 
in this context is how is the single chain dynamics affcx;tcd 
by chemical cross-links? It is our opinion the chemical 
cross-links should not be discernible from entanglements 
by just looking at the single chain dynamics. However, 
chemical cross-linking often leads to microphase syneresis 
[31], i.e., to separated microphases with the local mesh 
size being longer or shorter respectively than the average 
mesh size. The results derived below can also be applied 
to investigate such more complex situations [24]. 
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FIG. 1. The function 7(0,2/) of Eq. (2.14), which appears 
as a tinic-dopcndont prefactor in the approximate scaling law 
for k^nil{t) - k'niliQ) in Eq. (2.13). 

Let us further discuss Eq. (2.13) for the rather dilute 

case [xm = 0) without corrections from steric constraints 
first. Note that for very short times {t <C {k(.p)~^/^^^'^) 
the peaked function l{0,y) of Ec[. (2.14) shown in Fig. 1 
vanishes for all but the diagonal terms in the sum in 



Eq. (2.9). (Because we neglect end effects {kL ^ 1) 
and treat the polymer as homogeneous, we will speak of 
the diagonal term and write k^TZQ^{t) = ^(7fet)^/(0, 0) 
for k^'R.^^{t) .) As a consequence, the time decay of the 
structure factor is dominated by the diagonal element (s) 
of the sum in Eq. (2.9) and we have up to static terms 

^ ^ / P7eff2(i)\ crf(^fc7^(f(^)/2) 

,(.,t)ocV^exp( -fc^;,) - (2-15) 

This differs markedly from a simple exponential by its 
algebraically slow decay for long times but reduces to 

g{k, t) oc exp(-fc27^J-2(i)/6) (2.16) 

for the short times t ^ (A:£p)~''/'^7j7^; where it is sup- 
posed to be valid. For comparison, this function together 
with Eq. (2.15) and Eq. (2.17) is depicted in Fig. 2. We 
suspect that Eq. (2.16) may not be discernible in ex- 
periments, because it belongs to a time regime which 
falls within the crossover to the simple exponential ini- 
tial decay, i.e., the time {k£p)~^/^^^'^ ~ T/5-1 ~ Tj/Kk^ 
is of the same order of magnitude as the crossover time 
t* [33] to the initial decay regime discussed in Sec. Ill 
and in the Appendix. On the other hand, for suffi- 
ciently long times t » {kip)~^^^^ the function 1(0, y) 
may be replaced by the constant 7(0,0), because of the 
strongly oscillating static terms shown in Fig. 3. Note 
that in the limit of large kip "sufficiently long times" 
may be considerably shorter than the characteristic de- 
cay time 7^"^. Then, we can moreover neglect the static 
term 7^;^^(0) in Eq. (2.13). Hence, 7^;^^(^) is replaced 
by 7?,j|^(t) for all relevant times, and using the relation 

ir J2nm exp{ik^^TZim) oc (5(fcll), Eq. (2.8) reduces again to 
a simple exponential form 

.g(fc,i) cx exp(-fc27^^2(^)/4) (Hp > 1) . (2.17) 

The last result also could have been obtained immedi- 
ately from Eq. (2.8) for a purely transverse scattering ge- 
ometry with k = . This is what one expects, of course, 
since the longitudinal degree of freedom is suppressed and 
thus cannot contribute to inelastic scattering. Note that 
Eq. (2.17) is very similar to what one would have pre- 
dicted if the rigid constraint of constant contour length 
had been neglected altogether and isotropic motion of the 
local contour elements had been assumed. The end result 
in this approximation is obviously insensitive to the con- 
straint except that the mobility coefficient, which enters 
the calculation of TZQ^{t), is reduced as compared to the 
isotropic model. In other words, except for their differ- 
ent mobility an ensemble of randomly oriented scattering 
centers, which are constrained to move in two dimensions, 
cannot be distinguished from an ensemble of isotropically 
moving scattering centers by their dynamic scattering. 
This fact inspired us to use a simplified procedure to cal- 
culate the initial decay rate in Sec. III. It also explains 
the fortunate success of some of the "Gaussian" models 
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in fitting experimental data [19]. Let us emphasize, how- 
ever, that the reduction of the mobility coefficient due to 
the anisotropy of the undulations, which the "Gaussian" 
models fail to account for, is a crucial point that must not 
be neglected if one wants to obtain quantitative results 
and gain a real understanding of the underlying physics. 

Summarizing the above discussion we conclude that 
the time decay of the dynamic structure factor for a 
rather dilute solution {k^m ^ 1) of semiflexible polymer 
obeys 



g{k, t) oc exp ( — 



(2.18) 



(Here we have used /(0, 0) = r(l/4)/3.) For short times 
t if^^p)~^^^lk^ predict a deviation of the fimc- 
tional form of the decay from Eq. (2.18) according to 
g{k,t) cx exp(— 2r(i)/97r(7fci)3) and ultimately to the 
initial decay law derived in Sec. III. 
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FIG. 2. A comparison of the decay laws Eq. (2.15), 
Eq. (2.16), and Eq. (2.17) on a logarithmic time scale. In- 
dicated are the times U referred to in Figure 3. 

A stretched exponential decay as predicted by 
Eq. (2.18) has indeed been found in several experimental 
light scattering studies with semidilute actin solutions 
[23,6]. Assuming that the persistence length of actin 
is about 17 /im [34], wc conclude that k£p ~ 10^ and 
Eq. (2.18) should be well suited to describe dynamic light 
scattering from actin solutions. However, some fluores- 
cence microscopy studies [38] suggest that the flexibility 
of actin could be strongly length scale dependent and 
that actin may be much more flexible on the characteris- 
tic length scales probed in light scattering experiments. 
According to these results one would estimate k£p ~ 10"'^ 
and expect g{k, t) to cross over to Eq. (2.16) as discussed 
after Eq. (2.18). 

For the genuinely semidilute case, when fc^m — 1 and 
effects from steric constraints are not negligible, we have 
to consider the infrared cutoff > in Eq. (2.14). For 
IkTLm — {k^mY^^ of the order of 1 the structure factor 



decays in two steps as wc have mentioned after Eq. (2.14). 
We restrict the following analysis to the fast decay due 
to the single chain dynamics {t/TL^ = < !)• Then it 
is useful to split the integral in Eq. (2.14) into two parts, 
one of which runs from zero to infinity and the other 
from zero to Xm- A Taylor expansion of the integrand 
of the second integral renders fc^T^-g (f) to order 0{t^) in 
the simple form 
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FIG. 3. Static terms suppress contributions to the struc- 
ture factor which arise from correlations of distant contour 
elements. For several fixed times ti indicated in Figure 2, 
gnm{k,t) from Eq. (2.10) is shown as a function of kTZl,rn in 
the rod limit, i.e., for infinite persistence length £p. The os- 
cillations decay more slowly for larger times. For comparison 
the static structure factor of a rigid rod (solid line) has been 
included. 
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This expression allows one to extract the entanglement 
length L„i (and thus the mesh size ^„i) from experimental 
data and may eventually become useful in understand- 
ing the very interesting phenomena [31] that occur when 
semiflexible polymer solutions are gradually cross-linked. 
The same perturbation expansion has previously been 
claimed to account for hydrodynamic screening effects 
[32] . We do not agree with this interpretation: Hydrody- 
namic screening reduces the hydrodynamic correlations. 
Nevertheless, the mobility for the long modes certainly 
does not decrease to zero but saturates at the finite value 
2(±, as we have shown explicitly by use of Eq. (2.5). See 
also the discussion after Eq. (2.14). 



III. INITIAL SLOPE 

The results of the preceding section show that the 

structure factor g{k, t) of a solution of semiflexible poly- 
mers decays on a characteristic time scale 7^7^ by a 
stretched exponential law, Eq. (2.18). However, this ex- 
pression is not valid in the limit t ^ 0. The initial slope 
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of the structure factor does not become infinitely steep. 
Instead, as we explicitly demonstrate in the Appendix, 
g{k, t ^ 0) asymptotically approaches a simple exponen- 
tial law, exp(— 7^°''t) [33]. To calculate the initial decay 
rate 



(0) _ _ lng{k,t) 



dt 



(3.1) 



of the time decay of the dynamic structure factor we 
again suppose that Eq. (1.3) is fulfilled. A rigorous 
scheme for the calculation of the initial slope of the struc- 
ture factor for an ensemble of beads — coupled by a 
potential E{{rs]) and by the hydrodynamic interaction 
— may be found in the textbook of Doi and Edwards 
[9]. Here we assume that the stiffness of the molecule 
is not too large, so that the time scale separation be- 
tween solvent fluctuations and structural relaxation still 
holds, i.e., the characteristic relaxation rate of a bend- 
ing undulation of wavelength (^pA^)^/^ is supposed to 
be much smaller than the characteristic diffusion rate of 
solvent fluctuations of the same size. Then the hydro- 
dynamically driven small fluctuations about the equilib- 
rium conformation (where the potential forces vanish) 
are much faster than the structural relaxation of typical 
bending undulations and hence dominate the short time 
behavior of the dynamic structure factor. What is spe- 
cial about the case of a semiflexible polymer is again the 
rigid constraint imposed by the locally rodlike contour. 



If one would just follow the Doi-Edwards scheme naively, 
paying regard to the rodlikc structure of the molecule 
only through the static structure factor, the Oseen ten- 
sor would render each contour element two times as mo- 
bile along the contour as in the transverse directions, 
thus adding an artificial extra degree of freedom to the 
weakly bending rod problem, which is moreover weighted 
twice. On the level of a mere coimting of degrees of free- 
dom, one can thus say that — compared to the mobility 
matrix Eq. (2.2) — this amounts to overestimating the 
overall mobility of a contour clement by a factor of 2. If 
reversed, this argument tells that one may follow the sim- 
ple computional scheme of Ref. [9] if the result is divided 
by two in the end. This is indeed how we will proceed be- 
low. Of course, the procedure is not strictly rigorous: it 
docs not pay regard to the local anisotropy of the motion 
due to the rigid constraint. On the other hand, it takes 
the local reduction of the degrees of freedom properly 
into accoimt. From our experience with the problem of 
local anisotropy in the preceding section we may expect 
the result to be exact. 

Using the screened Oseen tensor in Fourier space. 



H{k) = 



1 



\ \k){k\\ 



(3.2) 



we arrive at the following integral for the initial decay 
rate 



(0) 

7fc 



l/a 



dq g{q) 



2r, J (27r)2 g{k) 



2^2 



■In 



q + k 



Mk^h 



(3.3) 



The static properties of the polymer enter the calcula- 
tion only via the static structure factor g{k) and through 
the factor 1/2. This implies a certain universality of the 
initial decay rate even in the semiflexible case. However, 
it turns out to be less universal than predicted by the 
Gaussian chain model. The main contributions to the 
integral come from large q. It actually diverges at the 
upper limit as a consequence of the short distance singu- 
larity of the Oseen tensor mentioned in the Introduction. 
This justifles the weakly bending rod approximation and 
for X <^ £p, L we may replace the static structure factor 
by its asymptotic form g{k) oc kr^. Again, we have in- 
troduced the lateral diameter a of the molecule to cut off 
the ultraviolet divergence. After doing the integral and 
dropping all but the leading order terms in ka one ob- 
tains in the dilute limit {k^h oo) the strikingly simple 
result 



7r 



keT 



k^ I ^ -\nka 



(3.4) 



Working with the full (rather lengthy) expression may 
be necessary in some applications, but the conceptual 
and practical significance of the result is best appreci- 
ated in the above limit. The deviation of Eq. (3.4) from 



the scaling law 7^°-' ~ k^ for a Gaussian chain is weak, 
hence it is useful to express Eq. (3.4) as a quasiscaling 
law 7^*^^ ^ k^^''^ with a wave-vector-dependent dynamical 
exponent z{k) given by 



z{k) = 3 



6 In fca — 3 
6 In fca — 5 



(3.5) 



This result is depicted in Fig. 4 together with a represen- 
tative curve for the more complicated semidilute case. It 
seems to account well for the deviations from the Gaus- 
sian chain prediction z = 3 seen in several experiments 
with polystyrene [35]. It is tempting to speculate that 
these are a signature of hydrodynamically enhanced lo- 
cal semiflexibility. It would be worth checking whether 
the initial decay of the structure factor measured by neu- 
tron scattering from solutions of polystyrene, DNA or any 
other intrinsically semiflexible polymer may be fitted by 
Eq. (3.4) with a reasonable value for a. Indeed, in the 
case of the biopolymer actin Eq. (3.4) already passed this 
test with remarkable success: It fits excellently available 
light scattering data with a = 5.4 nm [6], which com- 
pares very well to the value of 5.16 ± 0.3 nm obtained 
by much more elaborate methods [36,37] for two times 
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the transversal radius of gyration of the actin filament. 
Note that the resolution suggested by this comparison is 
far beyond the scattering wavelength, which was about 
an order of magnitude larger than a in the cited light 
scattering experiments. If the screening length has a fi- 
nite value, Eq. (3.4) and Eq. (3.5) have to be replaced by 
more complicated expressions, e.g.. 



(0) ksTk 
% = 



28- 30 ln{ka) + 6k-^^-^ 

- 6k-^^^^ avctan{k^h) - 30 fc"^^^^ arctan(fc^,,) 

- 3 k^a In (1 + fc-^C,-') - 15 In (l + fc-^e,"')] (3-6) 

for k^h ^ 1, but the qualitative structure is preserved. 
The main effect of screening is to flatten the increase of 
7^°'' and z{k) in the long wavelength limit (cf. Fig. 4). 
The above results for the dilute case can still be consid- 
ered a reasonable approximation for scmidilute solutions 
in a restricted range of scattering wavelengths X <^ ^m- 



dilute 

semidilute 




0.2 



FIG. 4. The effective (wave-voctor-depondont) dynamic ex- 
ponent z{k). The solid line is a plot of Eq. (3.5) for the di- 
lute case. The dashed line is a representative curve for a 
semidilute solution with screened hydrodynamic interactions 
(a = 200a). 



IV. CONCLUSIONS 

Dynamic scattering experiments allow for an accurate 
measurement of model parameters, if the number of these 
parameters is small and, of course, if the model is ad- 
equate within the range of wavelength and frequency 
probed. Scattering techniques are then sometimes more 
convenient than direct methods developed recently to in- 
vestigate the dynamics of individual polymers, such as 
flicker analysis [38] and microrheology [39] , which are still 
plagued by some technical complications [40] . They allow 
to probe in a quick and reliable experiment the internal 
dynamics of single polymers with the necessary statis- 
tical averaging inherently included. It was the purpose 



of our work to derive (within the common experimental 
accuracy) exact analytical expressions for the time decay 
of the dynamic structure factor for solutions of semiflexi- 
ble polymers. We succeeded in the case of Eq. (1.3), i.e., 
for scattering from the internal undulations of genuinely 
semiflexible polymers in solution. This is a generic situa- 
tion for neutron scattering from DNA or light scattering 
from actin solutions, just to mention two outstanding ap- 
plications. Our analytical results for the initial decay rate 
and the dynamic exponent for semiflexible polymers sug- 
gest that deviations of the dynamic exponent for (more 
or less) flexible polymers from its classical value z = 3 
are probably due to the local semiflexible structure of 
these molecules, and that the above analysis is therefore 
of some relevance also to scattering from rather flexible 
polymers. Moreover, these results lead to a convenient 
method to measure the microscopic lateral diameter a 
of a semiflexible polymer by use of the simple Eq. (3.4). 
From a fit of the result for g{k, t)/g{k, 0), 



exp 



iikt)i /r(i) ft 



1 



t 
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valid for kip » 1 and intermediate times t 



to experimental data one finds with good accuracy the 
bending modulus K (and the persistence length £p = 
n/ksT) of the molecules. For 1 the mesh size 

and the effective hydrodynamic screening length ^/j 
of the network can be determined. In the case of chem- 
ical cross-linking one has to be more careful since the 
development of microphases can cause the mesh size to 
vary considerably throughout the sample. For siiflicdently 
large kLm the terms containing tl^ can be neglected and 
the dynamic structure factor obeys Eq. (2.18). Note that 
the hydrodynamic screening length is closely related 
to the mesh size ^m, because the surrounding network 
disturbs the hydrodynamic autocorrelations of the indi- 
vidual polymers over distances longer than How- 
ever, one should keep in mind that — due to our ap- 
proximate treatment of the hydrodynamic interaction in 
Sec. II — in the effective friction coefficient used in 
Sec. II must be regarded as an effective parameter, which 
will in general differ from in its absolute value because 
it compensates for the deviation of the simple coefficient 
Eq. (2.2) from a more accurate description as Eq. (2.5). 
Finally we want to mention a hidden ambiguity of the 
above results. According to Eqs. (2.3), (2.4), and (2.12) 
a stiffening of the molecule (e.g., by a chemical manipu- 
lation) as well as an increase in concentration may cause 
a slowing down of the time decay of the structure fac- 
tor. An apparent increase in concentration may easily be 
caused by all kinds of syneresis in the course of chemical 
cross-linking, although the overall concentration remains 
constant. Therefore some caution is needed in the inter- 
pretation of experimental data. 
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For large p the parentheses are approximately equal for 
successive terms in the sum and wc can then replace the 
product of eigenfunctions by cos(p7r(s — <,)/L)/2. Sub- 
stituting the sum for t <C ti by an integral we arrive 
at 

^i-/(i)-^i'(o) = 



APPENDIX 

The eigenfunctions of the linear bending equation for 
the transverse undulations of a weakly bending rod with 
free ends at s = ±L/2 are given by [13] 

cos v^s/L cosh u^s/L 

u(p, s) oc + , ^ ' 

cosz/p/z coshi/p/2 

with Up w {2p + 1)tt/2 for p > odd and by the analoge 
expression with cos and cosh replaced by sin and sinh, 
respectively, for p even. Since we will only need the 
u{p, s) for large p ^ 1 and s L (far from the ends), we 
may drop the hyperbolic terms, which are small of order 
exp(— p7r/2), and write 



w(p,i)« (-) 



(p+l)/2 



cos 



for p odd and 



u{p,t)^{-Yl^sm(pn^ 



for p even. The transverse undulations and forces are 
expressed as 

ri{t) = 2 J2 rp{t)u{p, s), ft = jYl fp<P^ 



(We do not consider the zero mode.) The mode ampli- 
tudes decay as 



-t/Tp 



Here Tp := tl/p'^ and (vp)^ is determined by the 
equipartition theorem: 



2. ksTL^ 



k{-kpY 



After splitting Ti-if{t) into dynamic and static contribu- 



tions 



nj-^it) = 2 (r^(O)r^(O) - rj{t)r^iO)), 



insertion of the normal modes renders the dynamic term 



as 



u{p,s)u{p,(;). 



^/^cosM.-.)/Ll(l-e-'/'.). 



Introducing the scaling variables as discussed in the main 
text one now obtains Eq. (2.13). 

Finally we derive the initial decay within this ap- 
proach. In the limit t ^ 0, R^^{t) becomes small. In this 
case we may expand Eq. (2.8) to first order in R:j^^(t). 
Using Eq. (2.13) and Eq. (2.14) for large fc^„i and neglect- 
ing the static contribution (0) in the weakly bending 
rod limit, we get 



dx — I 

1 N ^ 



^ixka(n—m) 



l-(l-a:^)M!^! r^eos 



zka{n — to) 



After summation and neglecting kLx/2 against zcp :- 

zkL/2{t/TkLY^^ we have 

g{k,t Q) - g{k,0) = 



i/>,.^-l,^W,-f|(l-e-) 



SIX? Z(t) 



For short times t <^ TkL = Tfc-i (see also [33]) the z- 

integral reduces to 7r/2(f> and for large scattering vectors 
kL ^ 1 the static structure factor g{k,0) approaches 
n/ka, hence 

g{k » L-\ t^0)^^_ ^ksTj^s^ 



g{k,0) 



37ra 



The initial decay is simple exponential in time with the 
initial slope 

A more accurate expression for the initial decay rate 
is computed in Sec. III. The deviation is due to the 
approximate treatment of the hydrodynamic interaction 
in Section II. By comparison of the above result with 
Eq. (3.4) one can fix the phenomenological parameter 
for rather dilute systems (fc^m ^ 1), for which hydrody- 
namic screening is negligible. Then ^/^ becomes roughly 
equal to the scattering wavelength and no longer has the 
intuitive physical interpretation of a screening length. 



